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Stellar jitter from variable gravitational redshift: implications for 
RV confirmation of habitable exoplanets 
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ABSTRACT 

A variation of gravitational redshift, arising from stellar radius fluctuations, will introduce 
astrophysical noise into radial velocity measurements by shifting the centroid of the observed 
spectral lines. Shifting the centroid does not necessarily introduce line asymmetries. This is 
fundamentally different from other types of stellar jitter so far identified, which do result 
from line asymmetries. Furthermore, only a very small change in stellar radius, ^ 0.01%, 
is necessary to generate a gravitational redshift variation large enough to mask or mimic an 
Earth-twin. We explore possible mechanisms for stellar radius fluctuations in low-mass stars. 
Convective inhibition due to varying magnetic field strengths and the Wilson depression of 
starspots are both found to induce substantial gravitational redshift variations. Finally, we 
investigate a possible method for monitoring/correcting this newly identified potential source 
of jitter and comment on its impact for future exoplanet searches. 
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1 INTRODUCTION 

Space-based, photometric surveys have moved us into a new era of 
exoplanet discovery. The Kepler mission alone has identified 1235 
planetary cand idates, 54 of which lie within the habitable zone 
teorucki et allboi 1.) . However, radial velocity (RV) follow up is 
mandatory for planet confirmation. To do this for low-mass planets 
typically requires cm s~^ RV precision. 

However, astrophysical noise sources (or stellar jitter) due 
to spots, plages, granulation and stellar oscillations, for exam- 
ple, become an issue at the m s^^ level. These phenomena al- 
ter the shape of the stellar absorption lines, injecting spurious 
or systematic RV signals that may mask or mimic planetary sig- 
nals. There are numerous examples of this occurring in the liter- 
ature, e.g. HD 166435, HD 219542, TW Hya, LkCal9 and BD 
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l2008l ; lHuerta et alfcOOalFigueira et al.ll201 



^ Hy 



200ll:lDesidera et al.ll2004l:lHuelamo eFal] 



Spots and plages can induce RV signals of 1-100 m s~ 

jSaar & Donahue 1997). There have been many attempts to re- 
move this large-scale jitter using spot simulations, Keplerian mod- 
els, bisector a nalysis, harmonic decomposition, and Fourier anal- 
ysis etc. (see ISaar & Donahue! 1 19971 : iHatzesI |2002| : iDesort et all 
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I2OO7I; 'Dumusa ueetai] I2OII': 'Boisse et al.' '201 1'; 'Melo et 
2007; Bonfils e^lJl2007l :TBoisse et al. 2009; Oueloz et al. 200' 
iHatzes et all l2o1m Even with these removal techniques, active 
stars are still not ideal candidates for RV follow-up and are thus 
often left out of planet surveys. 

Unfortunately, even 'quiet' stars (those with little or no activ- 
ity) exhibit jitter due to granulation and stellar oscillations. Gran- 
ulation and stellar oscil lations can introduce typi cal RV signals of 
a m s~^ or higher (e.g. ISchriiver & Zwaadl200Ql . For this small- 
scale jitter, the currently implemented removal technique rests on 
adapting observational strategies to av erage out such noise (see 
Pepe et al.. 2011; Du musgue et am201lh . However, this technique 
is observationally intensive and does not provide information on 
the nature of jitter. 

In order to confirm low-mass planets around low-mass stars it 
is of utmost importance that we identify and understand stellar jit- 
ter, in all its forms. In this paper, we outline another possible source 
of stellar jitter due to variations in gravitational redshift. These vari- 
ations could potentially arise from fluctuations in the stellar radius 
or changes in the line formation height. 

In Section[2l we discuss the implications that variable gravita- 
tional redshift has on exoplanet confirmation/detection. We investi- 
gate measurements of solar radius fluctuations in Section|3] and in 
Section |4] potential mechanisms to produce stellar radius changes. 
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In Section [5] we discuss our findings and outline a possible solu- 
tion. 



2 GRAVITATIONAL REDSHIFT 

As photons escape from the deep gravitational potential well of 
their star they lose energy, causing a redshift in the obs erved spec- 
tral lines (e.g. Tavler 1994; Lindegre n & Dravinsll2003h . Any vari- 
ation in the gravitational redshift (GR) would result in a wholesale 
shift in the stellar absorption line positions which could mask or 
mimic a planetary signal. Unlike activity induced jitter, however, 
there are no resultant line asymmetries. Thus, jitter from GR vari- 
ations may prove more difficult to distinguish from the Doppler- 
reflex motion caused by a planet. 

Here, we estimate the magnitude of the potential RV signals 
introduced by such a GR variation. The observed velocity shift in 
the line centres, for a photon detected at an infinitely large distance 
from the stellar surface, is 

V (1) 

where G is the gravitational constant, M is the mass of the star, R is 
the radius of the star, and c is the speed of light. Given a fractional 
stellar radius change (at constant mass), 5R, the observed variation 
in velocity, SVgrav, would be 

SR 

SVgrav Vgrav- (2) 

For the case of the Sun, Vgrav is 636 m s"'^. If, however, the solar 
radius varied by only 0.01%, a distant observer would perceive a 
shift of ~ 6 cm s^^. The RV signal produced from the Earth around 
our Sun is only ~ 9 cm s~^. Thus a stellar radius change of ^ 
0.01%, by any means, could potentially mimic or mask a RV signal 
from an Earth-like planet. 

We are interested in the impact of GR jitter on all habitable ex- 
oplanet detections. The habitable zone scales with host star mass, 
as does the planetary RV signal. The RV signal from a habitable 
Earth-mass planet around spectral types FO-MO varies from about 3 
- 50 cm s~^. This means a much larger change in stellar radius is re- 
quired for lower-mass host stars to produce a SVgrav large enough 
to mask or mimic a planetary RV signal. This is shown in Figure[T] 
which displays the stellar radius change (expressed as a percentage) 
required to cause a SVgrav equal in magnitude to the RV signal due 
to a habitable Earth-mass planet orbiting stars of masses 0.3 - 1.8 

Mq. 

From Figure[T]we find, for M dwarfs, a relatively large stellar 
radius fluctuation (~ 0.1%) would be necessary to generate a sub- 
stantial enough GR variation to mask or mimic a habitable Earth- 
mass planet. Additionally, because a habitable planet orbiting a M 
dwarf has a period of days or weeks this radius change would need 
to occur over this timescale as well. For G2V and higher mass stars, 
we see a small (~ 0.01%) change in stellar radius, over a year or 
more, has the potential to induce GR jitter capable of masking or 
mimicking an Earth-like RV signal. Thus, GR variations are un- 
likely to pose a problem when planet hunting around M dwarfs 
due to the relatively large stellar radius change required over rela- 
tively short timescales. However, GR jitter could pose more signif- 
icant issues when searching around solar-like stars on account of 
the smaller GR variation required over far longer timescales. 

Because a radius change means that the stellar surface is re- 
ceding or approaching the observer, we have also examined the 
Doppler-shift imparted by this motion. Equating the amplitude of 




Mass (Mq) 

Figure 1. Smooth curve represents the stellar radius change (expressed as 
a percentage) required to induce a SVgrav equivalent to an Earth-twin RV 
signal. Circles represent, from right to left, spectral types: FO, F5, GO, G2, 
G5, KO, K5, and MO. Dashed curves represent stellar radius variations of 
50, 100 and 300 km. 

the RV signal due to a change in stellar radius (assumed to be a 
uniform contraction/expansion rate) over time to the corresponding 
Vgrav gives the characteristic timescale, tgrav, at which GR varia- 
tions begin to dominate, 

tgrav > (3) 

Thus, for low-mass stars we can rule out any radius fluctuations, 
such as stellar oscillations (e.g. the 5 minute solar oscillation), that 
occur on a timescale less than about 10 days as sources of GR jitter. 



3 THE SOLAR RADIUS 

Given that stellar radius changes can potentially mask or mimic an 
Earth-twin, we have searched for evidence of solar radius changes 
in the literature. Theoretically, there are a number of predictions of 
solar radius fluctuations ranging from ~ 0.00 1 to 0.0 1 % (where GR 
jitter bec omes an i s sue), l argely due to magnetic field interplay. For 
example, IStotherj l l2006l) argues using the virial theorem that en- 
ergy conservation requires a stellar ra dius decrease when magnetic 
activity increases. iFazel et ak Hugos') suggest that magnetic pres- 
sure interferes with convection and causes the size of solar gran- 
ules to decrease, lea ding to an apparent decrease in the solar radius. 
iMuUanetalJbOOTh also propose that magnetic pressure could con- 
tribute to hydrostatic support within the star, again leading to stellar 
radius fluctuations. 

Table [T] provides a summary of the most recent solar radius 
variation measurements that we are aware of (this table origi- 
nates from iThuiUier et J] jlOOSl) ; more recent results have been 
appended). While some authors report solar radius changes of 0.01- 
0.1%, the most recent space-based observations indicate that the so- 
lar radius is stable to within ~ 0.0001%. It is believed that the dis- 
crepancies between the reported solar radius changes may be due 
to atmospheric contamination of the ground-based observations. 
However, there is still an unexplained discrepancy bet ween results 
from stratospheric balloons and those made from space jSush et akl 
I2OIOI) . In conclusion, there is no concrete observational evidence 
of appreciable solar radius changes on the timescales comparable 
to the solar activity cycle. 
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4 POSSIBLE SOURCES OF GR VARIATIONS 

Given the lack of any clear observational evidence for solar radius 
changes of the magnitude required to raise substantial GR jitter, we 
now look at other possible mechanisms that may be important for 
other stars/spectral types. We note that each of these mechanisms 
will have their own jitter, however, it is beyond the scope of this 
letter to quantify these effects and we only consider the GR com- 
ponent. 



4.1 Magnetic Activity 

4.1.1 The Applegate Effect 

Quasi-periodic orbital period modulations, over a timespan of a 
few years, have been observed in several eclipsi ng binary sys- 
tems (e.g. ISaptista et al]|2003l ; iBorges et aHbOOSl and references 
therein); these ar e believed to be caused by the Applegate effect 
( lApDlegate|[l993) . According to the Applegate effect the stellar 
gravitational quadrupole moment changes due to the magnetic ac- 
tivity cycle, redistributing angular momentum within the star's inte- 
rior. The redistribution of angular momentum alters the oblateness 
of the star, thereby also causing stellar radius fluctuations. We at- 
tempt to estimate the magnitude of this stellar ra dius change, first in 
the close binaries from lLanza &Rodoncl ( fl99^ (excluding evolved 
stars) and then later extrapolate to exoplanet host stars. 

According to Rudiger et al. (2002), we can use the orbital 
period modulations, 5P/P, observed in binaries to estimate the 
change in gravitational quadrupole moment, SJ2, via 



(4) 



where a is the semi-major axis. The gravitat ional quadrupole mo- 
ment is alsorelatedtooblateness, / , and the stellar 

radius jSarnes & Fortnevll2003h via 



J2 



2 Re 



Re 



(5) 



where Re, and Rp are the stellar equatorial and polar radius, re- 
spectively. Differentiating equation |5] with respect to Req (where 
we assume Req ~ Rp) gives. 



SJ2 



25R 
3 R ■ 



(6) 



Thus, we can estimate the fractional stellar radius change from the 
observed orbital period variation. 



6R _ SP 



(7) 



Using equa ti on [7] and the binary star parameters from 
iLanza & Rodonq ([1999), we estimate stellar radius changes of 
0.001-0.1% for the rapidly rotating low-mass stars in these sys- 
tems. Unfortunately, such measurements of orbital period varia- 
tions do not exist for stars with rotation rates comparable to those 
of exoplanet host stars (which are generally slowly rotating). We 
therefore estimate the potential stellar radius change for these stars 
assuming that a fixed energy b udget is available to driv e the Apple- 
g ate mechanism (as done bv IWatson & MarshI 120101) . According 
to lWatson & MarslJ ( |2010|) . the change in energy, 5E, required for 
such orbital period modulations can be described as. 



1 G^MV 
6 Q.^R»Ms 



where Q. is the is the stellar rotational angular velocity and Ms is 
the mass o f the convective shell. 

As in IWatson & Marstil ( 1201 0^. the available energy budget is 
constrained to 10% of the stellar luminosity such that. 



5E < 



aLT 



(9) 



(8) 



where L is luminosity, T is the timescale of the gravitational 
quadrupole moment changes, and a is 0.1. Substituting equation 
|7]and|9]into equation[8] we find that the fractional change in stellar 
radius for a low-mass star over a year, due to the Applegate mech- 
anism, contributes less than a cm s~^ of GR jitter. Therefore GR 
jitter from the Applegate effect is most likely negligible in slowly 
rotating, low-mass exoplanet host stars. 



4.1.2 Inhibition of Convection 

The presen ce of a magnetic field has long been thought to inhibit 
convection ( lMakitjl958l : ISpruill2000l) . A decrease in the efficiency 
of convection is thought to result in the suppression of hot up-rising 
granules, possibly leading to an effective decrease in the observed 
stellar radius as the magnetic field is increased. Since the magnetic 
field over an activity cycle is variable, this mechanism may cause 
apparent stellar radius changes on the timescale of years. 

To further explore the impact that convective inhibition has 
on the effective stellar radius, we have simulated solar convection 
using the MURam code ( Vogler et al. 2005). The code solves large- 
eddy radiative three-dimensional magnetohydrodynamic equations 
on a Cartesian grid, and employs a fourth-order central difference 
scheme to calculate spatial derivatives. The numerical domain has 
a physical size of 12 x 12 Mm*^ in the horizontal direct ion and 
1.4 Mm in the vertical direction (see Sh elvag et al.ll201ll for de- 
tails). Continuum contribution functions were calculated in a nar- 
ro w (1 nm) band c entred at 630.2 nm using the method described 
in I Jess et al. l( l20TTI) . Since the total mass in the numerical domain 
is conserved and the gravitational acceleration is assumed constant 
throughout the domain, the change in height of the photospheric ra- 
diation formation in the model can be directly related to the change 
in the solar radius. 

Models of the solar photosphere were generated assuming an 
average magnetic field of 50, 100, 150 and 400 G. To reduce differ- 
ences introduced by stellar oscillations the contribution functions 
for each model were calculated at the same oscillation phase. We 
find that the 50 G model has a centre of gravity ~ 3 km higher 
than that of the 100 G, ~ 7 km higher than the 150 G, and ~ 40 
km higher than the 400 G model. This corresponds to a GR jitter 
ranging from ~ 0.3 - 4 cm s~^. Therefore, our simulations suggest 
that the inhibition of convection due to increasing magnetic field 
strength can be a possible source of GR jitter. 



4.2 Wilson Depression of Starspots 

It is well known that sunspots are depressed compared to the photo- 
sphere (the Wilson depression). The depth of the Wilson depression 
has been esti mated in many previous studies, ranging from 500 - 
2500 km (e.g |Suzuki|[l967ljBalthasar & Woehllll983l : lwatson et"ai] 
I2OO9I) . The spectral contribution of a sunspot with a Wilson depres- 
sion of 1400 km would be gravitationally red-shifted by ~ 1.3 m 
s^^ compared to the immaculate photosphere. However, as a spot 
contributes only a small amount of the total observed stellar flux, 
the net effect on the RVs will be much smaller. 

In order to estimate the amplitude of GR jitter from the Wilson 
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depression, we have created a model spotty star. Line-lists for spec- 
tral types G2V and K5V we re downloaded from the Vienna Atomi c 
Line Database (VALD-see pKupka et alJl2000l :lKuDka et al.l ll999h . 
Synthetic G2V and K5V spectra were then generated by convolv- 
ing these line-lists wi th an intrinsic line profil e with macroturbu- 
lence calculated from lValenti & Fischeil ( l2005h . wsini = km s^^, 
an instrumental broadening = 2.6 km s~^, and microturbulence = 
1.5 km s~^. The total synthetic spectrum for the spotty star was 
then modelled as the sum of the G2V spectrum (representing the 
immaculate photosphere) and the K5V spectrum (representing the 
cooler spots); we assumed a spot filling factor of 0.01. We tested 
two different depths for the Wilson depression, ~ 1400 km and 
3500 km; to do this the spectral contribution from the spots were 
shifted by ~ 1.3 and 3.2 m s~^, respectively (the shift expected due 
to OR). To measure the net RV shift induced by the OR compo- 
nent of the Wilson depression, we cross-correlated the non-shifted 
model spectrum with each of the two GR-shifted model spectrums. 
This process was repeated assuming a cooler spot, created using a 
MOV spectrum. 

For the model star with K5V spots, we find net RV shifts of 
^ 2 and 3 cm s^^ for Wilson depression depths of ~ 1400 km 
and 3500 km, respectively; for the model star with MOV spots we 
find net RV shifts of ~ 1 and 2 cm s^^. The model star with K5V 
spots produces larger net RV shifts because the spot component 
is brighter and therefore contributes more to the overall flux. The 
depth of the Wilson depression determines the magnitude of the RV 
shifts as a larger depth induces more GR jitter. Overall, for a quiet 
solar-type star, GR jitter from the Wilson depression is most likely 
a few cm s~^, which is not negligible when searching for habitable 
low-mass planets. We stress, however, that the presence of spots 
will introduce other sources of jitter. For example, since convec- 
tion is inhib ited within a spot, th e convective blueshift (typically 
lOO's m s" " MAsplundetal.i2000h is also inhibited. This would sig- 
nificantly shift the spotty spectrum relative to the immaculate pho- 
tosphere. 



5 DISCUSSION 

Changes in the observed spectral lines have the potential to mask or 
mimic the presence of a planet. In this paper, we have shown that 
a stellar radius fluctuation of even 0.01% can induce a GR varia- 
tion large enough to shift the centroid of spectral lines by ~ 6 cm 
s~^. We have also calculated a characteristic timescale, tgrav, over 
which stellar radius fluctuations must occur for GR jitter to domi- 
nant the RVs. From this we have found that GR jitter is a dominate 
noise source for stellar radius changes occurring on periods longer 
than ~ 10 days for low-mass stars. In light of this, we explored cur- 
rent solar radii measurements and potential mechanisms to produce 
such stellar radius changes. 

We find there is no universally accepted answer to the con- 
stancy of the solar radius. Most recent space-based measurements 
indicate a solar radius change of only 0.0001%. However, this is not 
in agreement with measurements from ground-based instruments, 
though there is reason to believe these may be systematically af- 
fected by the atmosphere. Missions such as PICARD and SDO may 
hold the key to high precision solar radius measurements in the fu- 
ture (The Picard Team et al. 2006; Pesnell et al. 201 1). 

We explored the origin of potential stellar radius changes. For 
slowly rotating, low-mass stars we find that the Applegate effect 
is most likely negligible. However, from our simulations we find 
changes in surface convection patterns, due to varying magnetic 



field strengths, alter the effective stellar radius sufficiently for ap- 
preciable GR jitter in the habitable planet regime. In addition, the 
Wilson depression of starspots was also found to contribute a sub- 
stantial GR shift when compared to Earth-like RV signals. There- 
fore GR jitter could be an important and hitherto unrecognised as- 
trophysical noise source. We note here that although significant, 
GR jitter may be overwhelmed by other line-altering mechanisms 
such as the Paschen-Back effect, which will be discussed in a forth- 
coming paper. 

Nonetheless, confronted with this new noise source, we be- 
lieve precise measurements of transit duration variations (TDVs) 
could be the key to removing GR jitter. As the stellar radius ex- 
pands or contracts, the transit duration would increase or decrease 
accordingly. Measurements of the TDV would allow us to calcu- 
late the change in stellar radius and thus monitor the GR variation 
at some level. 

This method does have its limitations as very precise photo- 
metric measurements are necessary to detect the potential TDV. For 
example, if the solar radius were to change by 60 km 0.01% R), 
then the transit duration of the Earth would differ by < 3 s. The 
lower limit of the TDV is set by transits with an impact factor of 
zero, while grazing planets would exhibit larger TDVs, and tran- 
sits may even disappear in extreme cases. Naturally, however, this 
would not enable the correction of GR jitter due to depressed spots. 

In conclusion, we have identified a new potential source of 
stellar jitter in the form of GR variations and outlined different 
mechanisms capable of generating cm s~^ variations. This may be 
important in the RV follow-up and confirmation of low-mass ter- 
restrial planets and Earth-like worlds. 
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6 Cegla et. al 



Table 1. Summary of observational measurements of solar radius fluctuations over an activity cycle, published primarily in the last ten years. Results are 
separated by location: G: ground-based, A: atmosphere-based, S: space-based. 
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